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SUMMARY 


This  report  presents  the  derivation  of  new  equations  for 
determining  the  changes  in  the  direction  of  a  ray  of  light  as  it 
passes  through  the  atmosphere  from  an  object  to  the  observer.  The 
equations  are  applicable  to  objects  inside  or  outside  the  atmosphere. 

Equations  are  also  derived  for  obtaining  the  topocentric 
distance  of  the  object  as  a  function  of  the  object's  height  and  the 
observed  zenith  distance. 


FORMULAS  FOR  COMPUTING  ATMOSPHERIC  REFRACTION 


FOR  OBJECTS  INSIDE  OR  OUTSIDE  THE  ATMOSPHERE 


I.  INTRODUCTION 

Existing  astronomic  refraction  equations  cannot  be  satisfac¬ 
torily  applied  to  objects  at  distances  up  to  a  few  thousands  of 
miles  from  the  earth  because  they  were  developed  for  a  special  use 
in  which  the  object  is  at  an  infinite  distance.  In  this  paper, 
equations  are  derived  that  are  applicable  to  objects  both  inside 
and  outside  the  atmosphere. 


II.  INVESTIGATION 

1.  Fundamental  Concepts.  Refer  to  Fig.  1  and  let 
A  =  observing  station 
S  =  position  of  the  object 
h  =  height  of  the  atmosphere  over  the  station 
C  =  earth's  center 

A  s  distance  from  the  satellite  to  the  station 
ZQ=  observed  zenith  distance 
AS=  curve  of  the  ray  path 
S.=  any  point  on  the  curve  AS 
e  =  astronomical  r. fraction 
AZ=  vertical  of  the  station  through  C 
R  =  refraction  to  be  determined 
x,y  =  rectangular  coordinate  system 
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The  observed  zenith  distance  ZQ  Is  the  apparent  zenith  distance 
of  the  object  when  the  ray  reaches  the  observer.  The  true  zenith 
distance  is  the  angle  ZQ  +  R,  *toich  the  ray  makes  with  the  vertical 
of  the  observer's  station  before  it  enters  the  atmosphere. 

The  direction  in  which  the  observer  sees  the  object  is  along 
the  tangent  to  the  curve  at  A.  The  origin  of  the  x  -  y  coordinate 
system  is  at  A,  the  y-axis  is  oriented  tangent  to  the  ray  path  at 
A,  and  the  x-axis  is  oriented  90°  to  the  increasing  zenith  distances. 


Let  N  be  the  point  where  the  normal  through  S  intersects  the 
y-axis.  The  desired  refraction  is  the  angle  NAS.  It  is  obtained 
by  equation 


sin  R  = 


(1) 


We  need  to  find  SN,  in  order  to  obtain  the  refraction  correction,  R. 

Because  here  we  consider  zenith  distances  less  than  79°,  R  is 
always  less  than  three  minutes  so  it  can  be  expressed  as  follows: 


R 


x 

A  sin  1" 


in  which 

x  =  SN  (2) 

We  find  x  from 

x  *  Je  dy  (3) 

so  that 


R  = 


1 

Asia  1" 


in  which  e  is  the  angle  between  the  y-axis  and  the  tangent  at  any 
point  of  the  ray  path. 

Assume  a  spherical  earth  with  the  atmosphere  arranged  in 
spherical  layers.  If  n  indicates  the  index  of  refraction  of  one 
layer,  n  +  dn  is  the  refraction  of  the  next  layer. 

If  in  the  first  layer  we  have  a  zenith  distance  Z,  the  next 
layer  is  Z  +  dc. 
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=ieg 


From  the  lav  of  refraction  we  have 


sin  (Z  +  de)  _  n 
sin  Z  n  +  dn 


(5) 


which  we  can  transform  as  follows 


so  we  obtain  for  de 


Then 


dn 

1  +  de  cot  Z  =  1  -  “?r 


dn 

de  =  -  tn  Z  q- 


f  dn 
e  =  -  Jtn  Z  -JJ- 

Insert  equation  (7)  into  equation  (4)  and  we  obtain: 

R  =  -  JI  tnZ^.dy 


(6) 

07) 


(8) 


solve  equation  (8)  we  need  an  expression  for  tn  Z  and  another  for 


To  develop  these  we  must  first  find  expressions  for  the  density 
of  the  atmosphere  (on  which  refraction  depends)  as  a  function  of  the 
height. 


2.  Density  of  the  Atmosphere  As  a  Function  of  Height.  The 
following  symbols  are  used: 

p ,  the  density  of  the  atmosphere  at  any  height 

p  f  the  density  at  the  earth's  surface  where  the  observer 
is  located  0 

h,  altitude  above  the  observer's  station 

n  .  index  of  refraction  at  the  observer's  station 
o 

n,  index  of  refraction  at  height  h. 

In  order  to  determine  the  refraction  it  is  necessary  to  have 
an  expression  for  the  density  of  the  air  as  a  function  of  the  height. 
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We  derive  an  empirical  lav  of  diminution  of  density  from 
observations  and  take  Into  consideration  the  fact  that  the  power  of 
reflecting  light  ceases  at  about  60  kilometers.* 

From  the  observed  values  of  the  density  at  different  heights, 
we  find  that  p  decreases  exponentially  with  altitude  following  the 
equation 


h 


in  which  h  is  a  constant, 
o 

The  values  of  density  p  up  to  20  kilometers  follow  equation 
(9)  with  accuracy.  For  heights  over  20  km  the  equation  is  less 
accurate,  but  still  sufficiently  accurate  because  the  refraction 
is  small  and  the  power  of  reflected  light  decreases  rapidly  at 
increasing  heights  (Fig.  2). 

The  constant  hp,  was  computed  by  weighting  the  observations 
proportionately  to  the  power  of  reflected  light  at  the  height 
of  observation.  The  value 


ho  =  9.2»*0  km 


and 


1  =  0.1082 

E  km 

o 


(10) 


Introduce  this  value  into  equation  (9)  and  we  obtain  for  the  density 
P 


„  -0.1082  h  (km) 

p  =  p  e 


(n) 


*The  basic  material  for  this  development  was  obtained  from 
The  Handbook  of  Geophysics  for  Air  Force  Designers,  Geophysics 
Research  Directorate,  Air  Force  Cambridge  Research  Center,  Air 
Research  and  Development  Command,  Cambridge,  Massachusetts  1957, 
and  from  available  balloon  observation  data. 
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From  Glandstone  and  Dale's  lav  the  values  of  the  Indexes  of 
refraction,  in  terns  of  the  atmospheric  density,  are  given  hy 
equation 


n  -  1 


P 


=  c  =  constant 


(12) 


in  which 


c  =  0.226 


Use  for  p  the  value  given  by  equation  (11),  and  the  index  of 
refraction  from  equation  (12)  can  be  expressed  as  follows: 


n  =  1  +  0.226  p 

o 


-0.1082  h 
e 


(13) 


in  which  the  height  h  must  be  expressed  in  km. 
From  equation  (13)  we  obtain 


,  -0.1082  h 

0.226  pQ  e 


1  +  0.226  p 

o 


-1082 

e 


0.1082  dh 


(110 


but  from  equation  (12) 


0.226  P  =  n  -  1 
o  o 


(15) 


and  because 

„  1 
0.1082  =  — 

ho 

We  can  rewrite  equation  (l4)  as  follovB: 


dn 

n 


-  * 

(nQ  -  1)  e  hQ  dh 
ho  Cl  +  (n0  "  1}  e‘ 


(16) 
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3.  Expression  for  ta  Z.  To  find  an  expression  for  tn  Z  ve 
must  consider: 

a.  The  maximum  height  of  the  stratus  over  which  the 
reflecting  power  can  be  assumed  zero. 

b.  The  maximum  zenith  distance  Z  • 

o 


For  (a)  we  find 


h  <  64  km 


and  for  (b)  we  adopt 


Z  *  75 
o 

With  these  values  for  h  and  ZQ,  the  maximum  value  of  Z  should  be: 

Z  =  Z0  -  5 


(17) 

(18) 

(19) 


in  which 


l  <  2?2 


We  can  then  express  tn  Z  as  follows: 

tn  Z  »  tn  ZQ  +  A 

in  which 

A  «=  f  (x) 

Esqpand  f(x)  using  the  McLaurin's  series.  Ihen,  we  have 


(20) 

(21) 


2  _3 

f(0)  +  x  f'(0)  +|  f*( 0)  +  g  f"(0)+..(22) 


Dilierentiating  equation  (20),  we  have 


i2  Z  >  « 
o 

After  the  second  differentiation 


(1  +  tn2  Z0)  g  »  f '(x) 


(23) 


^  (1  +  tn2  Z)  +  2(g)"  tn  Z  (1  +  tn2  Z) «  f*(x)  (24) 
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and  from  the  third  differentiation 


~  [l  +  tn2  zl  +  6  tn  Z  (1  +  tn2  Z)  ~  7^ 
,..3  L  J  hv2  dx 


dx 


dZ  ,3 


+  2  (g/  (1  +  tn2  Z)2  (1+3  tn*  Z)  *  f*(x) 


(25) 


To  find  the  — ,  ,  we  use  Snell's  law 

*  dx2  dx3 


nr  sin  Z  =  n  r  sin  Z.  =  constant 
0  0  ° 


Because  we  can  express  r  as  follows: 


*  -  r0  <1  ♦  f  ) 


(26) 


and  for  n 


n  =  1  +  (nQ  -  1)  e‘  HQ 


Snell's  lav  can  be  rewritten  in  the  form 


sin  Z 


_ n0  sin  Z0 


(27) 


Place 


we  obtain 


sin  Z  =  sin  Z.  (1  +  x) 


(28) 


(29) 
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It  follovB  from  equation  (26)  when  h  =  0 


r 


n 


n 


o 


then,  from  equation  (20) 


x  =  0 


and  from  equation  (29) 


Z  =  Z„ 


from  which  f'  (0)  =  0 


From  equation  (29)  we  obtain  for  x  =  0 

CsL  =  •tnZ° 

o 


(  ~  j  =  2  to  z0  +  tn3  ZQ 
dx2  o 


(  ~  )  =  -  3  tn  Z0  (1  +  tn2  ZQ)  (2  +  tn2  ZQ)  +  2  tn3  Zq 
v  dx3  t 

For  x  =  0,  equations  (23), (24),  and  (25)  become 


(30) 

(31) 

(32) 


f '(0)  =  (1  +  tn2  z0)  Qo 

2  2 

t*  (°)  =  C  “§  y  (1  +  tn2  Zo)  +  2  C  ^  ^  tn  Zo  (1  +  tn2  zc)  (33) 
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f"'(0)  .  [l  *  tn2  *o]  j_  (s§.}  +  6  <»  Zo(0y(i^ 
3 


d^Z 


♦  2  (  fjj  J  (1*3  tn2  Zc) 


Insert  the  values  of  equations  (30),  (31) ,  and  (32)  in  the  groups 
of  equations  (33)  and  ve  obtain  finally 


f '(0)  =  -tn  ZG  (1  +  tn2  ZQ) 

f"(0)  =  (2  tn  Z0  +  3  tn2  Z0)  (1  +  tn2  ZQ)  (34) 

f"'(0)  =  -  3  tn  ZQ  (1  +  tn2  Z0)  l2  +  tn2  ZQ  (7  +  5  tn2  Z0)  ] 


Nov,  introduce  the  values  of  equation  (34)  into  equation  (22)  and 
equation  (20)  becomes 


tn  Z  =  tn  ZQ  -  x  tn  Zq  (1  +  tn^  ZQ) 

2 

-  -  tn  Zft  (1  +  tn2  Z  )  i  2  +  tn2  Z_ 
2  o  o  L  o 


+  2-  tn  Z0  (1  +  tn2  ZQ) 
(7+5  tn2  ZQ)  ]  +  ... 


(2+3  tn2  ZQ) 
(35) 


If  ve  introduce 
A  =  tn  ZD  (1  +  tn2  ZQ) 

B  =  £  A  (2  +  3  tn2  Z^  (36) 

c  =  A  [l  +  i  tn2  Z0  (7  +  5  tn2  Z0)] 
equation  (35)  can  then  be  revrltten  as  follows: 

tnZ  *  tn  ZQ  •  x  A  +  x2  B  -  C  +  ...  (37) 
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This  equation  (37)  and  the  equation  (l6)  which  can  also  be  written 
in  the  fora 

4s  -  K  -  1)  e*  Eo  [l  -  -  1)  e-  EJ  to  (38) 

n  w  “o 

are  the  two  expressions  which  must  be  inserted  into  equation  #(8) . 
dy  remains  unknown  but  can  be  obtained  from 


dy  * 


dh 

cos  Z 


Its  value  is  found  later.  Let  us  first  solve  equation  (7) 


J 


dn 


k.  Solution  for  Jtn  Z  — .  Introduce  into  equation  (7)  the 
values  given  by  equations  (37)  and  (38),  and  ve  have 

e  =  -  tn  ZQ  -  (nD  -  1)  A  Jx  e"  E0  [l  -  (n0  -  1)  ^  ]  J” 

+  (no  -  1)  B  J  x2  e"  ^  [l  -  (nQ  -  1)  e_  hQ  ]  ^ 


(39) 


(nQ  -  1)  C  J  x3  e"  ^  [l  -  (no  -  1)  e‘  ho  ]  ^ 


The  solution  of  the  first  Integral  is 

-  in  32 .  r  _l_  t  -  °°  ~°2 

J”o  n  L  !  A  !n8,! 


m 


The  indexes  of  refraction  differ  frcm  unity  by  a  small  quantity 


n  <  1.0003 

Indicate  the  small  quantity  .0003,  by  y  and  we  have 

n  «  1  +  y  (*H) 
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in  which  y  has  the  values  given  by  equation  (26) 

_  h 

V  ■  (n0  -  !)  «  b0 

so  we  can  rewrite  equation  (40)  as  follows: 


- 

rx?'LV« 

2  n<j  n 


Neglect  terms  of  the 


J  2 

second  power  and  we  have 


m 


2  2 
po  - 

o  n  2  2 

2  n0  n 


=  Wo  -  v)  L1  ■  |  (Vo  +  v)  ] 


(43) 


v  is  the  value  for  h  =  0.  From  equation  (4l) 
o 


Y  =  n  -  1 
o  o 


m 


Introduce  the  values  of  y  given  by  equations  (44)  and  (42),  and 
equation  (43)  becomes 


n2  -  n2 
o 

_  2  2 
2  nQ  n 


h  r  -  -  1 

=  (nQ  -  1)  (1  -  e"  Ho)  [  1  .  |  (nQ  -  1)  (1  +  e  *o)  J  (45) 


With  this  value  the  first  integral  of  equation  (39)  is 

h 


-tn 


f  &  h 

z  Jll  -  (»o  -  1)  tn  ZD  (1  -  e  Hq)  .  1  (nQ  -  l)  (l  +  e"  r  )J  (**6) 


Solve  the  second  integral 


II 


J  x  e"  h0  [l  -  (nc  -  1)  e“  hj  ~ 


(*7) 
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in  'which  x  is  determined  hy  equation  (26) 


x  =  -  (n0  -  1)  + 


(nQ  -  1)  .  b 


no 


e  h«  +  — +  (nQ  -  1)  — e 


*o  no rQ  '“o 


noro 


Introduce  the  value  of  x  into  this  integral  (47)  and  we  have 


II 


=  -  (nQ  -  1)  J  e"  Eo  [i  -  (i^  -  l)  e"  Eo  ]  ^ 


(no  "  x)  J'  -  2  h. 


“o 

♦ 


e  ^  **0  Ll  -  (nQ  -  1)  e  Sq  ]  — 

Eo 


noro 


- E  r  - *  i  dh 

he  ho  Ll  -  (nQ  -  1)  e  Nj  J  “ 


(»»8) 


noro 


-  T  -  -  1  dh 

^  Ll  -  (no  -  1)  e  *o  J  “ 


h  e 


After  rearrangement  and  because  (nQ  -  1)  is  less  than  0.0003, 

n  +  (“o  +  ■  (n„  -  TJ  *  1  +  <"o  *  «  *  1 


and 


±  =  1  -  (n  -  1) 
n„  v  o  1 


we  obtain 


II  =  - 


(n0  -  1)  |’b  _  h  (n0  -  1)  |  h  -2  J 
— ~ —  J  o  e  hQ  dh  +  - J  o  e  *0  dh 
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1  i'h  -  k.  (go  -  jj  |h  -3  - 
+  n^rrJin  Jo  h  e  ho  ^  ■  iioho  J o  e  ho  dh 


(no  -  1) 


2  ,h  h 


j°h 


-3  - 

e  hQ  dh 


Integrate  to  obtain 


h 

hr 


-2  £ 


II  =  -  (n0  -  1)  LI  -  e  hD  ]  +  (no  “jj.  (i  .  e  ~  ho) 


ho  .  £  ,  h  'i  (n0  -  !)  [ 

*  Bl  -  «  hD  (  ^  .  1)  j  -  -TS—  U  -  e 


(no  -  l)2ho[  -3  £  .  ,  h  .  1 

- - - Ll  -  e  h  (l  +  3  r  )  J 

9  n  °  no 


vhich  is  the  solution  of  the  integral  of  equation  (47) 


Solve  the  third  integral  of  equation  (39) 


-  1) 


e 


-  *  i ^ 

ho\ 


III 


From  equation  (36),  for  the  maximum  zenith  distance  Z  =  75,  tn  Z 
3.7,  so  taking  tn  Z  =  4,  we  have  from  equation  (36) 


A  =  tn  ZQ  (1  +  tn2  Z0)  =  68 

B  =  |-  A  (2  +  3  tn2  ZQ)  =  1630 
but  because 

(nQ  -  1)  <  0.0003 


(50) 
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c f lty 


all  terms  containing  (no  -  ly  can  be  neglected,  and  we  have  for  the 
maximum  value 


B  (ne  -  l)3  <  0"003 


o  _  (no  -  l)h  h: 

x*  =  -2  — — — - + 


noro 


2  2 
^o  ro 


(51) 


Substitute  equation  (51)  into  (50)  and  we  obtain 


2(n0  .  l)  ,•  .  £  !  |  2  -  - 

III  =  -  - —  J  h  e  °o  dh  +  a  a  J  h  e  ho  dh 

“o  ro  ho  nS  r§  ho 

(“o  “  1>  f  2  -2  £ 

-  - X -  J  h  e  ho  dh 

n©  ro  ho 


(52) 


The  solution  of  these  integrals  is 


rh  .  h  2 

'Q  h  e  ho  dh  =  ho 


[x 


I* 


!0  h2  e"  ho  dh  =  h3 


f  »2  .■  ® «, -4- 1 


h  ,2  -i 

(53) 

»o  (£2  *  2  J  ♦  2)] 
hQ  **0 

(5h) 

"2  -  h2  h  1 

h0  (2  “2  +  2  £  +  1)J 

h^  ^ 

(55) 

Introduce  equations  (53) ,  (5*0  and  (55)  into  equation  (52),  and  we 
have 


2(no  -  1)  ^  T  -  £  ,  h  .1 
III - LI  -  e  hQ  (  —  +  l)  j 

“o  ro  ^ 
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(56) 


1  ^0  o  u  /  ** 

_  3-, 


,  h2  _  h  _x  ] 
(  “5  +  2  v  +  2)  j 


r.  -Sl,^2  .h  .vi 
°  ?  *>  1  - 


The  last  term  equation  of  (56)  has  no  influence  in  the  refraction 
correction.  By  introducing  this  term  into  equation  (39)  its  influence 
is  given  by 


,  B  ^°°  *  H  5°  .1  -  e’2  Ko  (2  *  3  S  ♦  1)  1  (57) 

3  "  »o  4  ^  *° 

The  maximum  value  corresponds  to  a  maximum  value  of  h  =  1^. 


It  follows  that 


f  _  2h  w2 


X-'  ho{2^-  -ho 


^  .  3  |  *  1)  ] 


max. 


=  O.lU 


For  tn  Z  =  6 


(Z  =  80?5) 
B  =  4218 


Retaining 


5o*o 


=  0.01 


(n.  -  1)  <  0.0003 


ve  obtain 


dR3  <  0*001 
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Similar  computations  for  the  second  term  give 


for  Z  =  75° 
for  Z  =  8l°5 


dR2  <  o'o6 
dB2  <  3*8 


(58) 


For  the  third  integral  we  make 


3  = 


x->  = 


h3 


1~3 
no  ro 


r  ,  -  £  r  .  h  -|  dh 

IV  =  J  x3  e  ho  LI  -  (no  -  1)  e  hj  ~ 


(59) 


IV  = 


h0  ^ 


r  3  -  £ 

JhJ  e  h 


dh 


(60) 


neglecting  the  second  term.  The  solution  is 


hJ  -  £  f 

IV  - - 2L  e  hQ  +  3  -2_3  L2  - 

*6*0  “0^0 


h  .2 


e  b°  ^  h2 


(61) 


Multiply  equation  (6l)  by  -C(np  -  1),  equation  (49)  ty  -A(np  -  1), 
the  first  two  terms  of  equation  (56)  by  B(n  -  1),  and  add  to  equa¬ 
tion  (46).  Then,  equation  (39)  can  be  written  as  follows: 


e  =  (hq  -  1)  tn  ZQ 
2  r 

+  (n0  -  1)  A  a 

(np  -  1)  A  hp 
*0*0 


(1  -  e  hp) 

Li  -  |  (nc 

-ll 

(np  -  l)2 

.  e  hp  J  . 

2 

r  * 

.  h  i 

.  h 

+  e  hp) 


-2  - 

A  (1  -  e  ho  ) 


(no  -  l)3 
3no 


[i 


] 
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ho  (no  -  l)3  r  1  1  r  -3  -  ,  hj 

*  — aU^oJ  Ll-e  h0(1O^)  J 


2  (Ho  -  l)2  ho  r  h 

-  B 

no  ro 


[l  -  e  ho  (Hq  +  1)  1 


(no  "  1)  h|  f  -  h  (  h2  h  >  1 

- 5—  B  L2  -  e  h  V  “5  +  2  C  +  2  J  J 


n§ 

(n0  -  1)  C  h|  /  5, 


o  ho 

o 


O  e  “0 


n 

“o  *o 


(62) 


3(no  -  1)  c  h|  r  -  r  r  h2  h  >)  *1 
- - - — — 2  12  -  e  ho  V  "2  +  2  h  +  2  J  J 


*0  ro 


h£  “  ho 


Collect  the  constant  terns  in  A,  and  we  have 


CA  =  (no  -  1) 


1"  ^  -  1  no  "  1  ho  ho  (^o  -  l)2  1 

l— s - ~ —  -  ~  ~  - J  (63) 


no 


noro  9 


Let  Cg  be  the  constant  terns  of  B.  We  then  find 


CB  '  2<“0  *  1)  L  ^-T  *  K  *  1) 


Crt  =  6 

C  n, 


nf  *§ 

(np  -  1)  hj 


“o  ro 


(64) 


o  r| 


_  h 

Let  Ei,  he  the  summation  of  terms  in  e  ^o;  Eq  of  the  e  ho,  and  E3 
_  3h 

be  of  the  e  ^o.  We  then  obtain 


2h 
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*1*  -  (“o  -  !)  z0  5^  - 


6  C  (n0  -  1)  h| 
“o  ro 


(65) 


+  2  B 


[(HO  '  1>‘ 


ho  (no  -  1)  hi 


ox  o  n|  i| 


n«r 


1 


("o  -1)  .  3 


A  H  2  (»o  -  !)  hn  Z* 


(66) 


E3  =  A 


hp  (pq  -  l)3 

9  no 


(i  +  l) 


(«7) 


n  -  — 

Let  be  the  sum  of  terms  of  ^  e  ^  and  M3  the  sum  of  terms  of 


h  -3  “ 

-  e  hjj.  we  then  have 
ho 


M1  =  (“o  *  1)  “TZ  A  +  2  B  [(no  -  l)2  +  (°°  g  **  —  ]  (68) 

Horo  ^  0  n£  rf 


“3-1  d*i> 


(69) 


and  letting  N  be 
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(70) 


3  K  -  1)  h3 


N  = 


»o 


r3 

*o 


C  - 


(qq  -  1)  ho  B 

— 


Substitute  equations  (63)  to  (70)  in  (62) 

«  =  (n,,  -  1)  te  Zp  [l  -  |  (Hq  -  1)]  A  CA  +  b  03  -  C  cc 

_  h  _2  h  _3  h 

+  E^  e  Hp  +  Eg  e  ^0  -  E3  e  hp 

h  ,  h 

„  h  “  V  „  h  "3  Z 

+  M1  hp  *  h°  ■  M3  hp  *  “° 

h2  _  h 


This  equation  (71)  gives  «  as  a  function  of  the  height  h,  vhich  we 
must  introduce  in  equation  (3)  to  find  x.  Tben 


o 


X 


{  (hp  -  1)  tn  Zp  [l  -  |  (np  -  1)  ]  +  A.Ca  +  B.Cb 


E2  J  e  \  dy  -  E3  J  e  hp  dy 


y  .3  h 


i7  V,  "  “  ry  .  .3  h 

+  “lJ0  Ee  6  *  “  M3Jo  i  *  *0  ^ 

2  _  & 

+  N  I  y  K  e  *9  ** 

Jo  h| 


(72) 
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The  first  integral  is 


J>- 

To  solve  the  other  integrals  we  need  an  expression  for  dy.  We  find 

^  =  cos  Z  (?3) 

To  find  -  -~~r  we  use  equations  (37)  and  (29). 
cos  £ 


tnZ  =  tnZQ-xA  +  Bx2-Cx3 
sin  Z  =  sin  ZQ  (l  +  x)  ^ 

By  division  we  obtain 


1 _  1  +  x  x  (1  +  x)  A  (1  +  x)  x2  B 

n  ~  '  -  . “ — -  +  V  ■ -  •• 


COS  Z  COS  Ze 


sin  Z(j 


sin  Z, 


h 
-m  — 


(7^) 


Because  in  all  the  integrals  the  exponential  e  (m  =  1,2,3) 

decreases  with  increasing  h  we  consider  only  terms  containing 
x^  and  neglect  the  other  terms  containing  x3  and  so  on. 

From  equations  (36)  we  obtain 


sin  Z„ 


B 


sin  Z. 


1  +  tn*  Zq 
cos 


4  (1  (2+3  tn®  Zj 


cos  Z, 


(75) 
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By  introducing  (75)  into  (7*0  we  have 


tn  Zq 


cos  Z  cos  Zn 


-  tntf  Zq  (1  +  tn2  Z0)  +  (76) 
cos  Z<,  2  cosz0  ^  0 


Multiply  this  equation  (76)  "by  dh,  and  equation  (73)  can  be  Witten 
as  follows: 


dy 


dh 

cos  Z0 


-*  <u>  +  2*a  tn2  80  (1  !  V  dh 

cos  “0  2  cos  Z0 


(77) 


Because  all  the  corrections  are  small  ones  and  decrease  with  in¬ 
creasing  height,  it  is  sufficiently  accurate  to  retain 


“e  ro 


-(=0  -  !) 


*2,  >4 


(78) 


Substitute  equation  (78)  into  (77)  and  we  have 


dy  = 


dh 


COS  Zq 


+  (Dq  -  1) 


tn2  Zq 

COS  Zq 


dh  - 


tn2 


_ ^0 

cos  Zo 


“oro 


dh 


+  |  tn2  Z0 


(1  +  tn2  Zq) 
COS  Zq 


(79) 


“o 


A  * 


Introduce  equation  (79)  into  each  integral  of  equation  (72)  and  we 
have  for  the  second  one 
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(80) 


i  £  •' t 


Ei  tn2  Z0  j*h  -  £ 

-  Z~F0  cos  Zo  J  o  h  *  ^  ^ 

,  Ei  tn2  Zq  (1  +  tax2  Zq)  fh  2  -  - 

*!  n|  rf  «.'i  Jo“  '  * 


dh 


Integrate  and  we  have 


El  f  N,  a?  »  [i*K-  »  ^  *»]  t1  ■  e  *•] 


+  El  t°2.Z°-j  [e~  i  (  s _  +  1)  -  l]  (8D 


noro  COB  *0 
2 


ho 


h  »2 


♦  [i D] 

nf  if  co®  Zq 


For  the  third  integral  we  have 


E<- 


Py  -  ^  [*  +  (°o  ~  1)  tn2  Zq]  jh  e"  ^ 

COS  Zq  •  o 


-y  .  gn 

j  e  h,,  dy  =  Eg 


dh 


_i_  ^  h  a, 

2  cos  Zq  r0  J  o 


(1  +  tn^)  ! 


-2- 


(82) 


♦  1  %  *D  ~  CO.  C"  WZ  ^  ^  “ 


the  solution  of  uhich  is 
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-2  ~ 


Eg  T  e  dy  =  Eg 


>*(«.-  J-)  h.  r  -2  ji  ,t 

coTZ^  ^  g”  L*e  K°+  XJ 


+  Er 


tn2  ^  if 


■2  r 


cos  ^  4  r0n§ 


[e  ^  (2  Eq  +  1)  -  l] 


(83) 


2  (1  +  tn2Zg)  jj3  r_  -2  -  h2 


*  g  E2  fac  — co-8^- 


2 

llorO 


[l  "  e  *>(2  h2  +  2  he+  ^ 


Solve  for  the  fourth  integral  and  we  have 


ry  -3  £  E, 

E  e  h0  dy  =-  — 1 
3  Jo  cos 


h-3- 

^  [1  +  (no  -  1)  tn2  Zq]  e  ho  dh 


E3  tn2  Zg  ph  _  3h 

+ - -r  e  h-j  .  h  dh 

n^  cos  Zg  J  o  ^ 


2  3  0  cos  Z0 


(8t) 


4  *i 


,  h 

u  h2  dh 


the  solution  of  which  is: 


'  e3  Z  ‘'3  «  -  §5®*  t1  *  <«»  -  1)  t»2  «.]  [•  "3  Eo  -  l] 

-2  „  .2 


9  n^  cos  Z0  L  ^  ho  J 

3E3  tn2  Zg  (1+tn2  Zg) 


>-3  [« -‘3  U  jf  ♦  «  ?,♦  «>] 


(85) 


54  ngrg  cos  Zg 
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For  the  fifth  integral  we  have 


py  h  -  -  _  [l  +  (dq  -  1)  tn2  Zq]  ph 

M1  -  e  h0  dy  =  K.  - J  *1 

1J°h0  0  ^  cos  Zo  Jo  h, 


Mi  tn2  Zo  ,h  h2  -  “ 

-  _i _  5  e  Oj  dh 

Ooro  cos  zo  J°  h° 


_  h 

e  hg  dh 


(86) 


+  3  &  Mi :  t»2  V  *  *3  5  ^ 

-  29  "  J°  ^ 


2  n2  rf  cos  Zq 


the  solution  of  which  is: 


^  Jl  Le  “  ■  jsri  l1  *  - 11  ‘“Ml1  -  •*  m^*  x>] 


.  iii  4  [a  .  e-  5.  (  $  ♦  2  -  2)1  (87) 

noro  cos  Zo  ^  L  h£  ^  J 


»  *£t  Z°(W  Z°’  4  [6  -  e-  O  £■  6  &  .  6 

2  nf  if  cos  Zq  **0  **0 


For  the  sixth  integral  we  have 


-  M3  Jro  ^  *’3  ^  ®  -  -  csh;  t1  *  (“c  - 1)  *\]  £  v'3  it  dh 
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”3 tng !»  rhh2e-3L 

n0r0  cos  Z0  J  0  h. 


dh 


(88) 


.  3  M3  tng  Z0(l  +  tn2  Zg)  fh  h3  .3  ^ 


2  nf  rf  cos  Z0 


f  h3  e  ^  dh 

Jo  \ 


the  solution  of  vhich  is: 


3h 


-m3  I*  |  e  he  dy  =  \  ^3  [l  +  (r^  -  1)  tn2  Z  1  Te"3  1^(3  |  +  l)-  ll 

J  vo  nc  9  cos  Zq  u  °J  L  ^  no  J 

M3  tn  Zq  .  r  .3  h  ^2  b  1 

-  ^r-o.-0  ^  -  «  >.  (9  lV  6  ^  «]  (89) 


3M3  tn\(l  *  r  .3  jl  h3  „a 

g_o - n|  L«  h»(2T  -3  ♦  27  2, 

162  rip  cos  Zq  0  hg 


+  18  S  + 


and  finally,  for  the  last  integral  we  have 


6)  -  6] 


;-y  h2 I  ....  _JL.  r,  .  /_  ,x  4._2n 


N  Jo  rs  e  ho  «y  -  sorz0  1.1  +  foo  - 1)  j0  .2 

“O 


If"  h! 

>J  j  o  V.2 


h 

e  h0  dh 


Ltn-23 .  f  h3  e-  i  dh 

n0r0  cos  Z0  Jo  ^ 


(90) 
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3  N  tn2  Z0(l  +  tn2  z0)  h  4  -  * 

- I  %  e  ho  dh 

J°  h2 


2  n?  rf  cos  Ze 


the  solution  of  vhich  is: 


ry  h2  -  - 

NJn^e  *0^  = 


h0  n  r 


o  h2  '  "w  COB  zo 


L^“  +  -  1)  tn  Zq 


,  ir  -*  h2  h  i 
2  ZjjJ  L2  -  e  ho  (  ^  +  2V  2)J 


N  tn2Z0  h2  r  -  ^  ,  h3  „  h2  ^  h  ^  r  "| 

+ - r  Le  he(  o+3  2  +  6  r  (91) 

n ere  cos  Zq  v,j  Vif  h- 


ho  bj 


3N  h|  tn2z^(l+tn2Z0) 


-  £  ^  .  h3 


|24  -  e  V  3  +  b  g 
2  Dq  Tq  cos  Zq  0  “0 


+  12  +  24  j*  +  24)] 

ho 


Use  the  abbreviation 


t  =  tn  Zp 


and  add  the  equations  (8l),  (83),  (85) ,  (87) ,  (89) >  (91)*  Then 

equation  (72)  becomes: 


x  = 


c  fly  =  {(^0  -  1)  t  [l  -  |  (no  -  1)]  +  A.CA  +  B.Cg  +  C.CC]  y 


+  E 


1 


bp 

cos 


(n,,-  l)t2Jl 


2  2 

V  bf  r  ' 

+  — - It 

n0r0cos  Zo 
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VO  1 1— f  u>|H  COM  rolK 


^(lrt2)  h|  r  -  £  ,  !  #  h  ,N  1 

3  73 - 7~  L1  "  e  h°(  2  +  h  +  x)  J 

nf  if  cos  Zp  Hq  no 


heLl  +  (tip-  l)t 
cos  Zp 


al 


(1  -  e”2  h0)+  i  E2  ^ ~ l_e  hp(2£ 
v  0/  k  d  n0r0cos  2pu  ^  n, 


r 


Eg  **<*  ^ *2>  j  r, .  „-2  s.  <2  a . 8  a  ,  D] 

2  co.  Z^n^rf  L  >«  “c  J 


Eq 


Jo-  [l  -  (n.  -  1)  t2  ]  [e‘3  h,  -  l] 


cos  Z, 


.2  2  3h  , 

E3  ,  110 _ rl  -  e  E;  (3  J  +i)l 


nor0  cos  Zo 


*tb  =3 ^v1  •  Jr [•*’ *“<9 5 +  6 1  *2)  - 2] 


*  %.  to^g;  [1  +  (».  - 1)  ‘2]  t1  - e  ^  ^  *  11  j  (92) 


%  t2  4_  [  -  ^  ( t?  |  +  2)  .  2] 

cob  Zp  nor0  L  0  v  ^2  h,,  J 

2  ^  nf  if  co.  Z„  L  *2  “I  "•  J 
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■f  coTj;  t1  *  K- 1)  <fl  u«'3  M3  ^  1)  -  l] 


m3  t2  h?  r  -3  r  h2  h  -) 
+  -  -1— r?-  L2  -  e  M9  j|  +  6  ^  +  2)] 


27  n0re  cos  Zg 


*  [l  *  K-  1)  t2]  L2  -  e'  5o(  j|  *  2  ^  *  2)] 


cos  Z, 


N  t2  hg  r  -  r  .  h3  '-2 


HgTg  COS  Zg 


r  -  v  ,  £  ,  £  *  h  ,1 

Le  i»  (  ^  *  3  ♦  6  -  *  6)  -  6\ 


♦  I  N  '±tlLLfl  [24  - h  *  * » s2  *  12  4  ♦  *  I  *  *1 
2  4  r|  cos  Z„  L  v3  ^  »l  *  J 


To  obtain  the  refraction  correction  we  divide  the  value  of  x 
by  the  distance  ts  ,  as  was  indicated  in  equation  (4) . 

In  order  to  consider  the  terms  necessary  for  Z  -  75°  we  must 
compute  the  values  of  E^,  Eq,  E^,  M^,  M^,  and  N.  We  find  for  Z  =  75° 

\  =  23  3" 

Eg  =  0.6 

e3  =  ° 

^  =  5.0 
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M3  =  1.0 


N  =  -0.05 

The  term  in  the  right  hand  of  equation  (65)  can  he  neglected  because 
its  maximum  value  reaches 


6  C  (hq  -  1) 


»? 

v? 


<  0^03 


and  each  coefficient  of  E-^  is  less  than  unity. 

Since  we  hold  as  maximum  zenith  distance  Z0  =  75° ,  R^fly  does 
not  reach  more  than  3  minutes  which  permits  us  to  retain 


A  =  y 


with  an  error  less  than  0.1  meter  for  the  maximum  value  of  A  (when 
the  object  is  at  the  limit  of  the  atmosphere) . 

Then,  the  final  formula  for  computing  the  refraction  for  an 
object  inside  the  atmosphere  is  as  follows: 


R  =  Aq  (n0  -  1)  tn  Zq  +  A1  (1^  -  1)  tn3  Zq  +  -  x)  tn5  ^ 


[(1  ♦  <«,-!)  tn2  Z0)  (X  -  e"  2„)] 


E!  ho  tn2  zo  « 

-  e 

HqXo  A  cos  Z„ 
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+ 1  ho  e2  rx 

2  t.  cos  Z0  L 


Mi  **<>  r 
A  cos  Z0  L1 


M3  * 

9  A  cos  Z0 


+  1) 


x] 


(93) 


+  N 


A  cos  Ze 


[l  +  K  -  x)  to2  zo]  u2  -  e'  ^(  75  +  2  r  +  2)] 


in  Tdiich 


Ao  =  +0.99827 
^  =  -0.00130 
Ag  =  +0.000006 

If  we  express  A  in  kilometers,  and  holding 

=  9.24  km 

h2 

no  =  0.0134  km 
*0*0 

1  h_  =  4.62  km 

2  ^ 

i  ho  =  1.027  km 
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The  height  of  the  station  hg,  does  not  affect  the  results  except  for 
extreme  heights.  If  this  is  so,  we  replace  h  and  r0  as  follows: 


where 


h  =  h8  -  ha 
Qoro  =  6372  ha 

hs  =  height  of  the  object 


ha  =  height  of  the  station 
both  hs  and  ha  are  expressed  in  kilometers. 

Equation  (93)  can  be  rewritten  as  follows: 


R  -  Ag  (n„  -  1)  tn  Zfc  +  Ax  (ng  -  l)  tn3  Z0  +  A2  (x^  -  l)  tn5  Zq 


9. 24  Ejl 
A  COS  Zg 


t- 


(no  -  1)  tn' 


?  *?  [1  -  .•  i] 


+  0.0134 


4.620  E 


tn  r  .  ^  ^  -1 

1  ^  [•  i.  ( *  ♦  «  - 1] 


A  cos  Zq 


bo 


2  r. 


A  cos  Zq 


1  -  e 


(<*) 


9.24o  (-  1  r  .  b  -) 

+  A  cos  Zo  ^  +  (“o  -  1)  to2  Zgj  Ll  -  e  bo  (  +  l)u 


1.026  M,  r  pir-3^h  -I 

+  r^"z~  -1  +  (”»  -  1>  ^oJ  L*  (3  -  +  1)  - 1] 
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9«24o  N  fi  + 
A  cos  Zq  l 


(no  -  1)  tn2 


®"  b®  (  ^  +  2  v  +  2)J 
*4  ho 


A©  =  +0.99827 


Ax  =  -0.00130 

Ag  =  +0.000006 


5.  Refraction  Viewed  from  an  Object  Inside  the  Atmosphere, 
The  refraction  viewed  from  the  object  is  the  angle 

a  =  NSA 


of  Figure  1.  We  see  that 


a  =  6  -  R 

The  angle  e  can  be  computed  by  using  equation  (71)  • 

Then  the  refraction  o  can  be  computed  as  follows: 


(95) 


-  £  9.2ko  e,  r  o  nr  ~  - 1 

=  %  e  ^ - i  J.  *  (no  -  l)  tn*  ZeJ  |_1  -  e  hJ 

A  COS  Zg 


,  E-,  tn2  Zq  r  -  -  ,  h  ,1 

+  °-013^  FcSTV  L1  '  6  ho  <  ^  +  «J 


-2  - 

+  E«  e  ho  - 


4.620 
A  cos  Zo 


n 

E2  [l  -  e'  ho] 


(96) 
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*  "i  K  e 


"  »0  . 


9.24o  M,  r. 


A  cos  Z, 


-  [l  +  (no  -  1)  tn2  Z0][l  -  e  h<,(  1) 


-  M- 


h  -3  h  1.026  Mo  T  o  T  ~3  ”  h  1 

-  e  ^ - 1  U1  +(n0-  1)  tn2  Zo][e  h,,  (3  £  +  1)  -  l] 

ho  ^  A  cos  Zj  JL  hg 


+  *  ^  e  "o-  /oofz,,  I-1  +  (n0‘  1>tn2  Z»I2  ■  e  ^  <  75  ♦  2  |L*  X 


As  always  /v  must  be  expressed  in  kilometers.  The  refraction  o  is 
always  less  than  R,  consequently 


e>R>| 

6.  Computation  of  the  Distance  £  of  an  object  inside  the 
Atmosphere .  When  the  distance  A  is  unknown  it  can  be  obtained  as 
follows: 


dy  = 


dh 

cos  Z 


From 


1  _  tn  Z 

cos  Z  sin  Z 

we  have  found 

1  =  1  .  x  <A  +  3  „2  tn2  Zp 

cos  Z  cos  ze  C0B  ^  5  cos  Zo 
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+  |  x\-QyV-  +  — —  (l  +  tn2  Z0)  + 

2  cob  Z0  cos  Zq  0 


•where 


r 


h  -m 


*>»  -  r  -=—  ^ 

^  “oro  ^ 


m  =  2,3  .... 

After  these  values  axe  replaced  and  Integrated  we  find 


*  2  „ 

tnZ^ 


_ h 

l 

h2 

tn2 

Zo 

/ 

cos 

Z0  5  norO 

cos 

Zo 

+  0 

- 

-  1)  ho 
“o 

tn2 

cos 

Zo 

Zo 

-  (°0  - 

.  3 

(no  -  1) 

h2 

tn2: 

Zc 

(1 

2 

Vo 

cos  : 

Zo 

+  3  ■ 

K  ■ *) 

4 

tn2 

Zo 

(1 

Vo 

COS 

Zo 

'  C08  Zq 

2  2 
\  ^  Zp 

i0r0  cos  Zq 
3Zo) 


2  T,2  ,.2 

no  ro 


tn2  Z0 
cos  Zq 


(1  +  -tsT  Zq) 


-  («o  -  i)  sr 


*o  tn2  Zo  h 


“e  cos  Z„ 


e 


(97) 
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e 


+  j  t?  K  -_L)  -  £  h  + © 

noro 


cos  Z„ 


(n©  ~  l)  2  ~  -  h  tn  Z0  o 

-  3  ^r”  ^  e  +  x)  cos  Zo  (!  +  tn  z0) 


7.  Working  Equations  to  Obtain  Coefficients  E1?  E,g  E^,  M^, 
M3,  and  N.  These  coefficients  can  be  computed  from  the  following 
working  equations : 


E1  =  -  (r^  -  1)  tn  Z©  +  [0^0680  -  0"00010  tn2  Z*  tn  Z©  (1  +  tn2  Z0) 


Eg  =  +  0"033  tn  Z0  +  0"008l  tn3  Z0 


E3  =  +  0"000005  tn  Z©  (1  +  tn2  Z©) 

(98) 

=  +  [o': 08^9  +  0"00044  tn2  Z©J  tn  Z©  (l  +  tn2  Z©) 


M3  =  +  0"0l68  tn  Z©  (1  +  tn2  Z©) 


N  =  -  0:000006  tn  Z©  (1  +  tn2  Z©)  (2+3  tn2  Z©) 

8.  Computation  of  (n©  -  l) .  The  refractive  Index  (n)  of 

standard  air  at  optical  frequencies  can  be  obtained  from  that  given 
by  Barrel  and  Sears 


rj  16.288  0.136 

K-  1)  107  =  2876.04  + - g—  +  — jp  (99) 

^  X 
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where 


X  =  the  light  group  wavelength  in  microns.  This  equation  agrees  to 
1  in  10®  over  the  visible  spectrum. 


(n<5  -  1) 


Pq  -  1 
1  +  at 


O.OOOOOOQ55  e 

1  +  Of  t 


where 


(100) 


=  refractive  index  under  ambient  conditions 

rig  =  refractive  index  in  dry  air  with  0.03  #  CO2  at  NTP  (0°  C,  760  mu 
Hg)  for  light  of  the  group  wavelength  employed,  as  calculated  here 

t  =  temperature  in  centigrade 

p  =  atmospheric  pressure  in  ran  Hg 

or  =  coefficient  of  expansion  of  air  (a  =  0.00 367) 

e  =  partial  vapor  pressure  in  mm  Hg 

Table  I  gives  the  vapor  pressure  corresponding  to  saturation  at 
various  temperatures 

Table  I.  Pressure  of  Saturated  Water  Vapor 


Temp. in  0°  C 

mm  of  Hg 

-5 

3.02 

0 

4.58 

5 

6.54 

10 

9.21 

15 

12.79 

20 

17.55 

25 

23.78 

30 

31.86 

35 

42.23 

1*0 

_ »•*> 
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Let  us  compute  an  example 


The  following  conditions  are 
observing  site: 

cp  =  latitude 

P  =  atmosphere  pressure 

t  =  temperature 

RH  =  relative  humidity  = 

\  =  effective  wave  length 

hs  =  height  of  the  object  = 

Zq  =  observed  zenith  distance  = 

ha  =  height  of  the  station  above 
sea  level 


chosen  as  being  in  effect  at  the 
38° 

760  mm  Hg 
-:-10°  C 

6o£ 

0.578 
13.96  tan 
70° 

0.1  tan 


Computation  of  the  Refractive  Index 


(n#  -  1)  107 


16.288  0.136 

2876.  OU  +  - 2 —  +  ■ -jp 


✓ 


(no  -  1) 


(na-  l) 

1  +  aft 


P 

760 


0.00000055  e 
1  +  cr  t 


\  =  0.578 
x2=  0.330625 
\4=  0.109313 


16.288 


0.136 

~f  , 

(n,-  1)  ID'  - 


. . .49.264 

1.244 

2876.04 

2926.55 


n,  -  1  =  0.000292655 
a  =  O.OO367 

From  Table  I,  for  60  $  R.  H.  and  temperature  =  10° 
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e  =  5.53 


so 


Hq  -  1  =  0.00028180 

(rig  -  1)"=  58"  125^ 

C omputation  of  the  Distance  £ 


h  =  hs  -  ha 


13.86  km 


*0 

B°  “  “ore 


O.OOIU50 


0.002176 


h_ 

ho 


1-5 


h 

cos  Zq 

1  Zq 

2  8  11  cos  Zq 


+ 


(no  -  1)  h 


tnfc  Zp 
cos  Zq 


-  (no  -  1)  \ 


2 

tn  Zp 
cos  Zq 


-  8o  •  ho  (^ 

-  |  s  h  (no 


1) 


2 

tn  Z0 
cos  Zo 


l) 


to2  Z0 
cos  Z0 


(l  +  tn2  Zo) 


1»0.5239 

-0.3328 

+0.0862 

-0.0575 

-0.0001 

-0.0024 
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+0 . 0022  tan 


tn  Z„ 


3  s»h"  ‘ 1(  (1  :  z»> 


1  o  2 

!'  2  h  8  (1  +  tn  Z°) 


+0.0062 


+  (n„  -  1)  ho 


2  h 

tn*  Z„ 


cos  Z. 


1  (n„  -  1)  h0sa  ?n  e  ho  (  r  +  1) 
cos  Zq  “0 

h  2 

-  3  (no  -  1)  h0s0  e  h0  (  £  +  1)  (1  +  tn2  Z0)  Z| 


cos  Zo 


+0.0128 


+0.0001 


S.  -0.0012 


A  =  40.2374  tan 


Computation  of  the  Coefficients  E1;  E2,  Mj_,  and  N 


Equation  (98) 


tn  Zo  =  2.747^77 
tn^Z0  =  7-5486 
tn^Zo  =  20.74 


2+3  tn^Z0  =  24-65 
tn  Zb(l  +  tn2^,)  =  23-49 


Ei  =  -(n0  -  1)  tn  Zp  +  [0"0680  -  0"00010  tn^l  tn  Z0  (l+tn^) 


Ex  =  -157-  H8 


E2  =  +0"033  tn  Zq  +  0'008l  tnJ  Zo 


E2  -  0"259 
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Mx  =  +[0"0l68  +  0" 00044  tn2  Z„]  tn  Z0  (l  +  tn2  Z„) 

Mj_  =  +  2"  072 


M3  =  +0*0168  tn  Zc  (1  +  tn2  Z0) 

M3  =  ;  0*395 


N  =  -0*00006  tn  (1  +  tn2  Z0)  (2  +  3  tn2  Z0) 


N  =  -0:035 

Computation  of  Refraction  R 
Equation  (94) 


A„  =  +0.99827 

cos  Z0 

=  0.34202 

A-l  =  -0.00130 

tn  Z0 

=  2.747477 

A2  =  +0.000006 

tn3  Z0 

=  20.74 

h  =  13.86 

tn5  Z0 

=  156.6 

h0  =  9-24 

h 

e‘  h0 

=  0.2231 

6  =  40.2374 

-3  - 

e  ho 

=  0.0111 

A  COS  Zg  =  13.7620 

l+Cno-  l)tn2z0 

=  1.0021 
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Ex  =  -157" 118 
E2  =  +0.259 


Mx  =  +2.072 
M3  =  +0.395 
N  =  -0.035 


h 


1 


e 


e 


2 


H 


-  e  h. 

=  0.7769 

ho 

=  -0.4422 

'3  h/h°  <3  *  <  D  -  1 

Ho 

_  h  .2  . 

-e 

«  -0.9390 

=  +0.3825 

(no  -  l)  Aq  tn  Zq 

159^422 

(ne  -  1)  Ax  tn3Z0 

-1.568 

(no  -  1)  Ag  tn?Z0 

+0.055 

9-24  e-l  r  x  P 

- -  _1  +  (n- -  1)  tn2 

a  cos  Z0  -  0 

1  r  -  -  “1 
^O-j  L^-  “  e 

-82.126 

2  h 

tn  Z.  1  ~  y. 

0.0134  Ex  - 2-  |_e  “0 

1  ^cos  Zc  L 

,  h  v  1 

<  V  11  •  ^ 

+0.511 

4.620  e2  >-  -  h  -| 

t.  COS  Zo  '-1  "  e  11(5 -1 

+O.O68 

43 


9.2k  mx  r 

A  cos  Zq 


[l  +  (n„-  1)  tn2  Z0j  [l  -  e  ^  (  j*  +  l)] 


1.026  Mo  r  .  r  h  h 

l1  *  ^  "  1)tn  ue  h0  (3  5^+  1)-  1^ 


A  cos  Zq 
9-24  N 


2  -  -!  r 


A  COS  Zq  L1  *  (no-  1)tn  zo_  L 


*>] 


+0.616 


+0.028 


+0.009 


77"  02 


R  =  77?  02 


Computation  of  the  Refraction  e 
Equation  (71) 

h 


-  =  1.5 

ho 

e"  ^  =  0.2231 

-2  h 

e  h0  =  0.050 

-3  - 

e  h0  0.011 

(«o-  1)  Ao  tn  Zq 

159-422 

:•  (n„-  1)  Ax  tn\ 

-  1.568 

+  (dq-  1)  A2  tn5Z0 

_  h 

+  O.O65 

+  V  »o 

-2  — 

-35-053 

*!*  Eg  6  ^0 

+  0.013 
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-3  - 
-  E3  e  h0 

0.000 

_  h 

+  Mj.  £  e"  *0 

+  0.693 

-3- 

-  Mo  -  e  ^ 

3*o 

-  0.007 

+  K  (  -  )2  e  4e 
*0 

-  0.018 

123?  537 

e  *  2'  03^54 

Computation  of  the  Refraction  a 

_  h 

®1  e  ho 

- 

35.053 

-  ■?—■■■■  [l  +  (no-  1)  tn2  zJ  [l  -  e"  b^l 

A  coa  Zo  4  J  «-  -» 

+ 

82.126 

E-,  tn2  2 L  r  _  h  .  -| 

+  1.0134  f - f- \l  -  e  h,  (  £  +  i)J 

A  cos  Zo  *>0 

- 

0.511 

+  Eg  e  l>o 

+ 

0.013 

4.620  lor  -  -  t 

-ACO.Z.L1-'  *] 

m 

0.068 

_  h 

+  M.  -  e  5>o 
^*0 

+ 

0.693 
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-0.6i6 


-  t  *  (B°- 11  to2 <lt-  •" *  <  i* «] 

A  cos  Zq  “o 


_3h 

MX  i  •  *• 


-0.007 


1.026  M 


-3  - 

2  [i  +  (no-  1)  tn2  zj  [e  ^  (3  |  +  1)-  l]  -0.028 


A  cos  Zp 


B  -* 


+  H  (  —  )  e 
*0 


-0.018 


9-240 
A  cos  Z0 


[l  +  (np-  1)  tn2  Z0]  [2  -  e  *0  (  H2  +  2  j|  +  l)]  -0. 

bp  00 


009 


46252 


o  =  46*.  52 
a  =  c  -  R 
a  =  46*52 

Astronomical  Refraction 


e_  =  (np  -  l)  A0  tn  Zp  +  (np-  l)  Ax  tn3  Zp+(i»p-  l)A2tn5Zp 


«a  *  157*91 


The  set  of  values  obtained  are  as  follow: 


«a  =  2'37*91 
«  «  2^3*54 
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R  =  1  17-02 


CT  =  46.52 


i'i7"02  >  l'oi"8 


9.  Refraction  Viewed  from  an  Object  Outside  the  Atmosphere. 

The  foregoing  equations  are  correctly  applied  to  an  object  inside 
the  atmosphere  for  zenith  distances  less  than  75°  and  with  the 
height  of  the  atmosphere  of  65  kilometers,  which  may  be  regarded 
as  the  limit  beyond  which  the  air  does  not  produce  any  appreciable 
refraction.  For  such  conditions  we  have  expanded  tn  Z  in  series. 

An  object  outside  the  atmosphere  can  be  treated  as  follows: 

In  Figure  3  let  S  be  the  object.  There  is  a  point  S0  at 
which  tiie  ray  coming  from  S  reaches  the  upper  layer  of  the  atmosphere 
after  which  following  the  curve  S0A,  it  reaches  the  point  A. 

The  x  and  y  coordinates  have  the  same  meaning  as  before. 

At  the  limit  of  the  atmosphere,  formula  (94)  gives  R0 . 

«a  -  THS  is  only  the  astronomic  refraction. 

The  refraction  we  are  looking  for  is  Rg  =  TAS. 


From  Figure  3  it  follows  that 


sin  Rg 


sin  Rg 


fo 

»0 


(101) 
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Fig.  3*  Fundaaent&l  concepts  In  relation  to  Fig.  1. 
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NT  tn  e» 


but 


euad 


*s  - 
NT  =  NP  +  PT 


80 


then 


Xg  =  Xo  +  PT  tn  «a 


sin  R-  =  ^2.  +  —  tn  «a 

^  A,  A.  * 


(102) 


Because  Be  or  Rq  is  a 


angle  of  about  3  minutes  for  Ze  =  750 


so 


W  =  As-  A0 


*b  As“  Ao 

sin  Re  -  —  + _ tn  a, 

*s  Aa 


(103) 


The  refraction  Rq  can  be  obtained  from  equation  (94)  without  considering 

_  h  _3  h 

the  terms  which  contain  e  v  ®nd  e  Eq,  because  at  the  height  of 
8e,  they  become  extremely  small.  Since  we  can  write 


*a  *  (n©-  1)  tn  [l  -  |  (ug-  1)]  +  A.CA  +  B.Cg  +  C. 


Cc 


We  can  write  equation  (93)  as  follows: 

*.  ■  *•  *  T  tsX  (1  +  (B»-  x>  **  *•  -  >] 
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but 


and 


Xg  *  MT  tn  «a 
NT  =  NP  +  PT 


80 


then 


Xg  =  Xg  +  PT  tn  «a 


sin  Rg  «  +  —  tn  «• 

*  *S  As 


(102) 


Because  Rg  or  Rg  is  a  snail  angle  of  about  3  minutes  for  Z0  =  75° 


w-  *g-  *o 


so 


A.- 


Xo 

sin  Rg  =  —  + 

*■  A. 


tn  «. 


(103) 


Itoe  refraction  Rg  can  be  obtained  from  equation  (94)  without  considering 

.  h  _3  h 

the  terms  which  contain  e  hg  ®nd  e  Eg,  because  at  the  height  of 
80,  they  become  extremely  small.  Since  we  can  write 

*a  -  (flg-  1)  tn  [l  -  |  (ng-  1)]  +  A.Ca  +  B.Cjj  +  C.Cfc 


We  can  write  equation  (93)  as  follows: 

"o  ■  -  *  T  EsSSt.  (1  +  11  *2  •»  ■  ^7*  0 
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(u*) 


♦  ifs  “p  .h  ». 

2  COB  z0  £  COB  Z0  AQ  9  cos  Z0 

Because 


*o  =  Ag  •  Bln  Rq  =  A  • 


It  follovB  that 


*•  -  A  *  h  L  f1  *  K-  «  ^  z„.  h?~;Zt’)] 


(105) 


+  i  Eo  ,^£  +  M, 

2  d  cos  Zq  * 


M- 


ho 


cos  Z0  ^9  cos  Z0 


Indicate  by  Z  all  the  terms  of  the  right  hand  free  of  A  and  ve 
have 

Xo  =  «*Ae+  2  (106) 

Nov,  introduce  this  value  Xq  into  equation  (103)  sod  ve  obtain 


®s  =  *a 


«a 


vhich  simplifies  to 


Rs  - 


(107) 


(108) 


because 


e 


e 


a 
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Equation  (108)  shows  that  the  refraction  outside  the  atmosphere  is 
obtained  by  introducing  a  correction  to  the  astronomical  refraction, 

E 

by  an  amount  given  by -  . 

As 

Since  can  be  any  distance,  for  an  object  at  infinity 

=  0 

As 

consequently, 


*e  =  €a 

which  is  the  astronomical  refraction. 

10.  Computation  of  the  Distance  ^  of  an  Object  Outside  the 
Atmosphere.  Equation  (97)  is  not  good  now  for  computing  the  distance 
of  an  object  outside  the  atmosphere  because  it  was  developed  to  be 
used  inside  the  atmosphere.  For  a  height  over  65  kilometers  the 
following  equation  must  be  used: 


where 


A 


ro 


1  +  cos  6 

cos  £  e  COS  Oy  -  <j) 


(109) 


e  *  Zo  +  «a  -  z 

ea  »  astronomical  refraction 

XL. 

Z  *  arc  sin  (^j-  sin  Zq) 

«a  “  (°o“  1)  &o  tn  Ze  +  (x^-  l) A1tn3Z0+  (x^ -l) A2tn5Z0 
Y  -  i  (Zo  +  ea  +  Z) 
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a  alleys  is  a  small  angle  (a  fev  seconds)  so  we  can  assume  a  ■  0  for 
computing  the  distance  A  ,  vith  sufficient  accuracy  for  obtaining  R 
or  a.  To  obtain  an  error  dR  ■  0"01  the  error  6A  in  A  reaches  the 
following  values  for  a  zenith  distance  Zq  «  70°: 

H  =  100  km  is  6A  =  200  meters 

H  =  1000  km  is  6A  *  10,000  meters 

If  it  is  required  to  obtain  a  vith  higher  accuracy,  in  equation  (109) 
the  value  of  a  computed  from  equation  (ill)  must  be  introduced. 

Then,  the  refraction  R  can  be  obtained  by  using  the  following 
equation: 


R  =  (np-  l)  Ap  tn  Zp  +  (np-  l)  Aj  tn3  Zp 


+  (np-  1)  Ag  tn5  Zp 


9-2*  %  T  o  1 

+  A^TZp  L1  +  (bp  -  1)  tn2  Zp! 

B,  tn2  Z0  4.62  E„ 

-  0.0134  — - -  +  - - 

A  cos  Zp  A  COS  Zp 


+ 


[1. 

A  cos  Zp 


1.026 


r 

m3  Li  + 


(np  -  l)  tn2  Zp] 

(np-  1)  tn^Zp] 
Acoa  Zp 


+ 


18.46  R 
A  cos  Zp 


(110) 


11.  Computation  of  the  Refraction  Angle  q  of  an  Object  1 
the  Atmosphere.  This  refraction  is  the  angle  a.  It  can  be  oh 
by  using  the  following  equation: 


outside 


Trrrr. 


[l  +  (np-  1)  tn2  Zp]  +  0.0134  ^ 

A  cos  Zp  L  ^  A  cos  Zp 
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4.62  t  Mi  Ll  +  (n©  -  l)  tn  V 

_ i  -  9-24 - 


£  COS  Zq  a  cos  Zo 

r  2  1 

Mo  Ll  +  (i^  -  l)  tn  Z0J  l8. 1*6  N 

+  1.026  — — — “  -  A  AA0  7 
A  COS  Zq  A  COS  Zq 


(111) 


Example  of  Computation 

We  assume  the  same  data  as  used  for  computing  R  and  a  for  an 
object  inside  the  atmosphere,  and  change  h8  to: 


as  before 


and  assuming 


the  constants  used  are: 


h  =  100  kilometers 

8 

hg  =  1,000  kilometers 
Zg  =  70° 
ha  =  0 


■l- 

-  157"  H8 

E2  - 

+  0.259 

“l  = 

+  2.072 

of 

II 

+  0.395 

N  = 

-  0.035 

\>  = 

6370.06  km 
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Computation  of  the  Distance 

xu 

sin  Z  =  __2 -  sin  Z* 

(l  +  s)  ^ 


H  =  100  ton 

H  =  1,000  tan 

s 

O.OI56984 

0.1569840 

sin  Ziq 

0.939692 6 

no 

1.0002818 

sin  Z 

0.9254296 

0.8124204 

Z 

670  43'  59// 6 

54°  19'  59-0 

Zo 

70° 

70° 

«a 

2'  37-9 

2  37-9 

y 

68  53  18.7 

62  11  18.4 

e 

2  18  38.3 

15  42  38.9 

i  e 

l  09  19.2 

7  51  19.4 

1  +  cos  e 

1.999  1868 

1.962  6407 

cos  -y 

O.36018  36 

0.466  5650 

A 

277.579 

2123.206 

A  .COS  Zq 

94.938 

Y  -  a 

58°  53*03*6 

62°  Ul'i6*3 

cos  (y  -  a) 

0.360  2519 

0.466  5741 

A 

277*5264 

2123.1698 

A  cos  Z0 

94.920 

726.167 
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Computation  of  the  Refraction  a 

H  g  100 _ H  =  1,000 


Ex  [l  +  (no-  l)tn2  zj 

-  9.24 

A  cos  Zo 

+15-422 

+2.016 

Ex  tn2  Zj, 

+  0.0134 

A  cos  Zo 

-O.167 

-0.022 

4.62  Eg 

-0.013 

-0.002 

,A  cos  Z0 

9-24  Mx  [l  +  (no-  1)  tn2  Z0] 

-0.202 

-0.026 

A  cos  Zq 

r- 

M3  Li  +  (no-  1)  tn2  ZpJ 

+  1.026 

A  COB  Zq 

+0.043 

+0.001 

18.48  N 

A  COS  Zo 

+0.006 

+0.001 

a  = 

+15" 083 

I.968 

Now,  introduce  these  values  of  a,  and  the  new  values  of  A  ere: 

H  =  100  km: A  =  277-5264  km 

H  =  1000  km: A =  2123.1698  km 

These  new  values  do  not  alter  the  preceding  value. 

Then,  the  results  axe  as  follow: 

H  -  100  km  H  =  1000  km 
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